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Potential Condensed Fuel for the Milky Way 

M. E. Putman^ 

ABSTRACT 

Potential condensed clouds of gas in the Galactic halo are examined in the 
context of the recent models of cooling, fragmenting clouds building up the bary- 
onic mass of the Galaxy. 582 high-velocity clouds (HVCs) are defined as the 
potential infalling, condensed clouds and the sample's spatial and velocity distri- 
bution are presented. With the majority of the hydrogen in the clouds ionized 
(~ 85%), the clouds at a distribution of distances within 150 kpc, and their 
individual total masses below 10^ M0, the total mass in potentially condensed 
clouds is 1.1 — 1.4 X 10^ M©. If the tighter distance constraint of < 60 kpc is 
adopted this mass range drops to 4.5 — 6.1 x 10^ M0. The implications on the 
condensing cloud models, as well as feedback and additional accretion methods, 
are discussed. 

Subject headings: Galaxy: halo — intergalactic medium — galaxies: formation — 
cooling flows — Galaxy: formation 

1. Introduction 

The range of stellar ages and metallicities in galaxies like the Milky Way indicate fresh 
star formation fuel must fall in throughout their lives (e.g., Rocha-Pinto et al. 2000; Renda 
et al. 2005). The gas accretion process has traditionally been thought to proceed via shock- 
heated halo gas from the intergalactic medium cooling and falling in to feed the star formation 
process (e.g.. White & Rees 1978; White & Frenk 1991). Recently these models have been 
extended from all of the gas within a certain radius collapsing monolithically, to including 
fragmentation as the hot gas cools, forming pressure supported warm clouds (e.g.. Mailer 
& Bullock 2004, hereafter MB04; Kaufmann et al. 2005; Sommer-Larsen 2006). Models 
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which include fragmentation do not have the "over-coohng problem" the monolithic collapse 
models have. In another words, all of the gas does not cool and fall in at once, a large 
fraction remains in the halo in a warm/hot phase and excessive feedback is not necessary to 
explain the observed baryonic mass of the galaxies. MB04 predict at any given time during 
a Milky Way-like galaxy's recent evolution, several thousand condensed clouds with a total 
mass on the order of 2 x 10^° can be found within a ~150 kpc radius. These clouds are 
pressure confined by the hot gaseous halo remaining and should currently be found in the 
Galactic halo as evidence of this ongoing gas accretion. 

Likely candidates for these condensed infalling clouds are the high- velocity clouds (HVCs) 
of neutral hydrogen surrounding our Galaxy. Oort (1966) was the first to propose this type 
of origin for HVGs and MB04 discuss the similarities between the HVGs and the condensed 
clouds in their model. HVGs range in size from ultra-compact (< 20 arcmin^; e.g., Briins & 
Westmeier 2004) to hugely extended (1800 deg^; Wakker & van Woerden 1997) and have typ- 
ical peak column densities of approximately 10^^ cm~^ (Putman et al. 2002). Their velocities 
generally extend from 90 km s~^< |VL5i?| < 450 km s~^, or —300 km s~^< Vqsr < 300 km 
s~^. Recent progress on the distances to HVCs allows origins such as the condensing cloud 
model to be seriously considered and constrained. The direct distance constraints involve 
looking for absorption lines in the spectra of halo stars at known distance and generally 
provide lower limits on the order of > 5 kpc, but also include upper limits for 3 clouds of 
< 10 kpc (summarized by Wakker 2001; Thom et al. 2006). Indirect distance constraints 
include deep HI observations of systems similar to the Milky Way and Local Group (< 160 
kpc; e.g., Pisano et al. 2004; Zwaan 2001), Ha observations indicating the HVGs are being 
ionized by photons escaping from the Milky Way (< 40 kpc for some clouds; e.g., Putman 
et al. 2003a), and constraints on the properties of the compact HVCs (CHVCs; 6 < 1°) 
when subject to the extragalactic ionizing background radiation (< 200 kpc; e.g., Maloney 
& Putman 2003). In addition, recent work surveying the environment of M31 does not find 
the M31 analogs of CHVCs beyond 60 kpc (Westmeier, Briins & Kerp 2005). All of these 
distance constraints place the clouds within the extended Galactic halo and appear to be 
consistent with the distances expected for the condensing, infalling clouds in the models. 

Given the recent developments both theoretically and observationally regarding the 
gaseous distribution about galaxies, the time is right to address the properties of potential 
condensed clouds currently in the Galactic halo. This paper addresses the observational 
constraints on the HVG population and how this can be put together with the models 
to form a consistent picture of condensed, infaUing clouds feeding the Milky Way's star 
formation. The selection of high-velocity clouds is presented in the next section, followed 
by the properties of this sample in the context of the condensed cloud model. Finally, the 
results are discussed and summarized. 



-3- 



2. Data: HVC Selection 

Potential condensed clouds were selected from an updated version of the all-sky HVC 
catalog of Wakker & van Woerden (1991; hereafter WvW91). The catalog has been up- 
dated by including clouds from the catalog of HVC components by Morras et al. (2000) for 
declinations < —23° and includes clouds with |Vlsr| > 90 km s~^ (see Wakker 2004 for 
more information on the catalog properties). The updated WvW91 catalog of 626 clouds is 
not exactly high resolution (0.5° at best), but covers the entire sky with a detection limit 
of 2 — 3 X 10^* cm~^ (Av = 25 km s~^). The only selection criteria for the HVCs to be 
potentially condensed halo clouds was the exclusion of clouds associated with the Magellanic 
System (e.g., the Magellanic Stream and Leading Arm; Putman et al. 2003b) and the Outer 
Arm Complex. The Outer Arm Complex is a large low latitude structure that is consistent 
with being a warped section of the outer Galactic Disk, or a high-z spiral arm (e.g., Wakker 
& van Woerden 1997). The remaining 582 clouds have a mean velocity relative to the Galac- 
tic Standard of Rest of Vgsr — -43 km s~^. The clouds in the condensing cloud models are 
predicted to have a range of velocities as they form and move through the halo, but should 
have a net infall in agreement with this mean Ygsr (Bullock pers. comm.). The clouds have 
a total HI flux of 993,557 Jy km s~^. The largest contributors to this total flux are Complex 
C (209,590 Jy km s'^) and then Complex H (98,040 Jy km s"^). 



3. Results 

The spatial distribution of the potential infaUing HVCs is shown in Figure 1 with the 
symbol representative of a positive (star) or negative (triangle) Vgsr cloud. The clouds have 
Vgsr's between approximately +/— 300 km s~^. The clouds are distributed throughout 
the sky with positive and negative velocity clouds intermingled. The largest concentra- 
tions of clouds are towards the anti-center in the southern Galactic hemisphere and around 
/ = 260° in the northern Galactic hemisphere. These clouds are thus found looking away 
from the majority of Galactic disk, however there is also a number of clouds found in both 
hemispheres from I — — 45°. The clouds with Vqsr more negative than -100 km s~^ are 
also concentrated between I — — 45° and around I — 180°. 

The total mass of the entire population of potentially infaUing clouds as a function of 
average distance is shown in Figure 2. This plot assumes the fraction of the cloud that is 
detectable as neutral hydrogen is 10.5% of the total mass of the cloud; 70% of the cloud is 
hydrogen and 15% of that hydrogen has cooled to be detectable as HI. This 15% is justified 
by assuming the extragalactic ionizing radiation field is the dominant factor in ionizing the 
clouds (e.g., Maloney & Putman 2003; see Section 4.3). The HI mass was calculated using. 
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Mhi = 2.36 X 10^ I (Jy km s-^) D2(Mpc), with I = 993,557 Jy km s"^ for all of the clouds 
selected according to the previous section. The total mass, Mtot, is then Mij7/0.105. With 
this constraint. Figure 2 shows that if the high velocity HI flux is at an average distance of 
80 kpc the total cloud mass is ~ 2 x 10^° Mq, at 40 kpc the total mass would be ~ 6 x 10^ 
Mq, and ~ 9 X 10^ Mq at an average distance of 20 kpc. 

There are several clouds that contribute significantly to the total HI flux in halo clouds of 
993,557 Jy km s^^, and it is more realistic to place the clouds at a distribution of distances 
rather than an average distance. For instance one of the larger clouds which is part of 
Complex A has a flux of 52,210 Jy km s~^ and is known to be between 4.0 — 9.9 kpc (van 
Woerden et al. 1999). In addition, total individual cloud masses between 10^^^ Mq are 
considered the most likely initial mass for each cloud in the models given the constraints on 
the formation of the clouds (i.e., the ability to fragment and the conduction limit), cloud 
survival (i.e., Kelvin- Helmholtz Instability), and cloud motion (i.e., ambient drag) (MB05; 
Kaufmann ct al. 2005). Placing the clouds at a random distribution of distances < 150 kpc 
and ensuring that the total mass of each individual cloud remains below 10^ Mq (again with 
the HI 10.5% of the total mass) results in a range of total masses in condensed clouds between 

1.1 — 1.3 X 10^ Mq. All of the range of masses are presented at the 95% confidence level and 
were found after running over 100 random distributions of distances. The lower limit on the 
mass of inidvidual clouds is left open, as the clouds are expected to be disrupted as they 
come into close proximity with the Galaxy. An example of the distribution of distances and 
total individual cloud masses generated is shown in Figure 3. If we adopt the tighter distance 
constraint of < 60 kpc for all clouds, this range of total mass decreases to 4.5 — 5.7 x 10^ 
Mq. If all of the clouds were actually given the range of masses between 10^~^ Mq and 
kept within 150 kpc, the total mass in potentially condensed clouds goes up somewhat to 

1.2 — 1.4 X 10^ Mq. If all clouds are given the same mass there is no direct correlation 
between the clouds' resulting distances and their observed GSR velocities. This might be 
expected if the small clouds represent distant clouds not yet affected by the Galaxy's halo 
medium or gravitational pull. 

Clouds that have not yet been detected by existing surveys will have low HI fluxes and 
most likely small masses, but if numerous they could significantly change the total mass 
in condensed clouds in the halo. The effect these yet undetected clouds may have on the 
total mass has been tested by extrapolating the HVC flux distribution function of log N(S) 
= —1.44 log S + 3.91 (Wakker 2004), where N(S) is the number of clouds with a given flux 
S, measured in 10 Jy km s~^ bins. This distribution holds to approximately 25 Jy km 
s~^ for the updated WvW91 catalog and the fact that it does not continue further is at 
least partially due to the completeness limits of the catalog, as also found in the southern 
HIPASS HVC catalog (Putman et al. 2002). When the HVC flux distribution function is 
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extrapolated to continue to 1 Jy km and these low flux clouds are assigned with the 

same random range of distances within 150 kpc, the total mass in clouds increases only 
slightly to 1.1 — 1.4 x 10^ Mq. Keeping the clouds within 60 kpc and adopting this flux 
distribution function results in a total mass range of 4.5 — 6.1 x 10^ Mq. Thus, with a flxed 
ionized component, the potentially existing lower flux clouds will not add significantly to 
the total mass in condensed clouds. If a model is adopted in which the smallest flux clouds 
(< 100 Jy km s~^; a flux cut that encompasses the majority of the compact HVCs; Putman 
et al. 2002) are at closer distances (< 20 kpc) and the rest of the HVC population extends 
to 150 kpc, the range of total mass in condensed clouds is 0.9 — 1.1 x 10^ M©. If on the 
other hand the clouds with small HI fluxes (< 100 Jy km s~^ again) are presumed to be at 
distances greater than 60 kpc (but < 150 kpc) the total mass in condensed clouds increases 
to 1.3 — 1.6 X 10^ Mq. Finally, if all of these clouds were actually given a mass range between 
10^~'' Mq (somewhat unrealistic, as it places the clouds with fluxes of 1-2 Jy km s^^ between 
150-200 kpc), the total mass reaches a similar 1.4 — 1.6 x 10^ Mq. Keeping all of the clouds 
within 150 kpc requires lowering the bottom of the cloud mass range to 5 x 10^ Mq and the 
total mass in condensed clouds lowers to 1.3 — 1.5 x 10^ Mq. 



4. Discussion 



The total mass of potentially condensed clouds in the Galactic halo is limited to be at 
most approximately 1.6 x 10^ Mq in all of the above cloud distributions that are consistent 
with the distance constraints and keeping 15% of the hydrogen in each cloud neutral. This 
upper mass limit is set by extrapolating the flux distribution function to encompass yet 
undetected HI clouds, placing the clouds at a random distribution of distances below 150 
kpc, and constraining their individual total mass to be below 10^ Mq. This mass is over a 
factor of 2 lower if the clouds are kept within 60 kpc. The 3 main factors affecting this total 
mass are the distances, the limit on the total mass of individual clouds, and the percentage 
of the cloud that is neutral. Each of these factors are discussed here, followed by a discussion 
of the impact of these results on the gas accretion models. 



4.1. Distances to HVCs 

As discussed in the introduction, there are now a large number of constraints on the 
distances to HVCs. All of the constraints are consistent with placing the clouds within 150 
kpc and some of the constraints place the clouds at closer distances. If the HVCs extend 
only out to 60 kpc and not 150 kpc and the clouds are kept at 15% ionized, the total mass 
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in condensed clouds changes from 1.1 — 1.4 x 10^ M© to 4.5 — 6.1 x 10^ Mq. If we refer back 
to Figure 2, this type of total mass would place the majority of the HI flux in the halo at an 
average distance of approximately 21 kpc or 12 kpc, respectively. In either case, the upper 
limit on the total mass in condensed clouds is dependent on how many clouds are at largest 
distances. This will be constrained further with ongoing distance determination programs 
(e.g., Thom et al. 2006). The limit on the total mass in clouds if they are within 60 kpc is 
consistent with the findings for M31 of 3 — 4 x 10^ Mq (assuming 10.5% HI again) in clouds 
within 25 kpc of this galaxy (Thilker et al. 2004). Closer distances for the clouds in the 
models may infer the cooling times are longer than initially presumed and/or the densities 
of the fragments in the outer Galactic halo are not high enough. 



4.2. Individual Cloud Masses 



The mass of individual clouds as the condense within the hot halo depends on a number 
of factors and is not currently tightly constrained. MB04 find a range of masses are suitable, 
with 10^~^ Mq being the most likely given the constraints of conduction, evaporation, Kelvin- 
Helmhotz instability, and ram pressure drag. Kaufmann et al. (2005) and Sommer-Larsen 
(2006) find a similar range of masses for individual clouds. Tidal disruption is a factor that 
could disrupt the clouds within approximately 13 kpc and lead to some small clouds that no 
longer have typical masses in this range (MB04). 10^ Mq is therefore adopted as the upper 
limit on the total mass of the individual clouds, and since some HVCs may lie at distances 
below 13 kpc and represent condensed clouds that have been disrupted, the lower limit on 
the mass of individual clouds is left open, as indicated in Fig. 3. If the small clouds, or 
lower flux clouds, are presumed to lie at distances below 20 kpc (rather than extending out 
to 150 kpc) while the rest of the clouds are allowed to extend out to 150 kpc (as long as 
their mass remains below 10^ Mq) the total mass m clouds drops to 0.9 - 1.1 x 10^ Mq. 
One could also argue that the small clouds should be more distant however; simply based 
on their angular size. Placing the small clouds at distances greater than 60 kpc, but within 
150 kpc, increases the total mass in condensed clouds to 1.3 — 1.6 x 10^. The middle ground 
of placing the clouds at a range of distances appears to be the best approach. In contrast to 
the small clouds, the largest HVC complexes are unlikely to be beyond 20 kpc, and this is 
kept the case by constraining a cloud's total mass to be below 10^ Mq. 
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4.3. Neutral Fraction 

This paper uses the nominal value of 15% of the hydrogen in each cloud being neutral 
based on the impact of the ionizing flux from the extragalactic background light (10^ photons 
cm^^ s~^; Maloney & Putman 2003) and the majority of the clouds lying at distances between 
60-100 kpc. The possibility that a larger percentage of most clouds are ionized should be 
considered. A larger ionized component could be due to many of the clouds lying at closer 
distances and being subject to the ionizing radiation from our Galaxy (e.g., Putman et al. 
2003a) or via coUisional ionization as the HVCs interact with the remaining hot halo medium 
(e.g., Sembach et al. 2003). It is also possible some of the small clouds are more distant and 
a lower fraction of the cloud has cooled to be observable in HI. If only 1% of the hydrogen in 
each cloud is detectable in HI, the total mass of clouds in the halo would reach approximately 
1.7 - 2.1 X 10^° Mq with the clouds within 150 kpc and 6.8 - 9.2 x 10^ M© with the clouds 
within 60 kpc. A 1% neutral component is possible for some clouds given the discovery of 
highly ionized HVCs (Sembach et al. 1999; Fox et al. 2005), but is unlikely for the majority 
of the HVCs given the results of pointed Ha measurements (e.g., Tufte, Reynolds & Haffner 
1998; Putman et al. 2003a). In fact the current hmited Ha measurements generally show a 
larger fraction of neutral material than ionized for HI HVCs. There are some HI HVCs that 
show evidence for extended ionized components (e.g., Sembach et al. 2003; Haffner 2005), 
but others do not (Putman et al. 2003a). Though proximity to the Galaxy may result in 
more of a cloud being ionized, this may also be offset by the clouds closer to the Galaxy 
harboring higher density material. Estimating the actual fraction of ionized material relative 
to neutral is difficult given the limited pencil beam sightlines used to probe the ionized 
component. 15% neutral is a reasonable estimate based on the theoretical predictions and 
current observational constraints. Future sensitive Fabry-Perot Ha observations should help 
to clarify the full extent of the ionized component of HVCs. 

The fraction of each cloud that is neutral will also depend on the amount of metals in 
the gas. As with the direct distances, there are very few HVC metallicity determinations. 
The metallicity determinations for the clouds in the sample presented here are almost solely 
limited to the giant Complex C and generally range from 0.1 - 0.3 solar (e.g., Collins et al. 
2002; Tripp et al. 2003). The model of MB04 uses a metallicity of 0.1 solar, so consistent 
with, but on the low end of the limited HVC metallicity estimates. If they included more 
metal-rich gas, the gas would cool more efficiently at lower densities and larger radii from 
the Galaxy. The HVCs would then be expected to lie at even larger distances (~ 200 kpc) 
which seems unlikely given the HVC distance constraints. Higher metallicities in the MB04 
model would also most likely result in an increased total mass of condensed halo clouds. 
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4.4. Implications on Gas Accretion 

Disk galaxy formation models generally predict the existence of a hot halo with a bary- 
onic mass of a few x 10^° M0 (e.g., MB04; Fukugita & Peebles 2005; Sommer-Larsen 2006). 
Hot halo gas has recently been detected in the vicinity of our Galaxy (Sembach et al. 2003; 
Rasmussen et al. 2002) and around other spirals (Pedersen et al. 2006) in possible support 
of these models. This hot gas gradually cools as clouds and fuels the galaxy's star formation, 
but the specifics of the process vary by model. The model of MB04 predicts that cloud 
formation and infall balances at early times, with the balance of condensed clouds in the 
halo approximately 2 x 10^° Mq. The analysis of the observed HI halo clouds presented 
here indicates the total mass in potentially condensed clouds in the halo is at least an order 
of magnitude below this. The < 6 x 10* M0 range found here for the halo clouds within 
60 kpc is consistent with the simulations of Kaufmann et al. (2005) and Sommer-Larsen 
(2006). The reason for less mass in condensed halo clouds than found by MB04 may be due 
to the clouds falling in rapidly after they are formed, and thus less clouds are visible in the 
halo at a given time. If clouds fall in rapidly, feedback from the Galaxy may need to be 
considered to keep most of the halo material in a warm-hot phase and not over-produce the 
baryonic mass of the Galaxy. A Galactic fountain is one possible feedback mechanism in 
which the hot gas from multiple supernovae rises into the halo (e.g., Bregman 1980). The 
Galactic fountain could inject a large amount of enriched hot gas into the lower Galactic 
halo which may then mix with the massive hot halo material and cool as the intermediate 
velocity clouds (IVCs) found around our Galaxy. IVCs are much closer to disk (0.5 - 2 kpc) 
than HVCs and are also of higher metalhcity (0.5-2 solar; Wakker 2001). 

The limited mass in condensed clouds forming and falling on to the Galaxy indicates 
additional accretion methods are necessary to explain the Galaxy's observed baryonic mass. 
Several of these additional accretion methods are directly evident as our Galaxy destroys 
smaller galactic systems such as the Sagittarius dwarf and Magellanic Clouds. The Magel- 
lanic System itself will eventually bring on the order of 10^ Mq of HI into the Milky Way 
(Putman et al. 2003b). The Magellanic System is an example of a less frequent, large accre- 
tion event and could also potentially disrupt the process of cloud formation at the current 
epoch. Finally, though cold accretion is unlikely to dominate at low redshifts and for galaxies 
as massive as the Milky Way, this method of gas accretion in which the incoming gas is never 
heated to the virial temperature of the galaxy halo, may also need to be considered (e.g., 
Keres et al. 2005). In any case, in the model of Keres et al. the smaller galaxies currently 
being accreted by the Galaxy obtained the bulk of their mass via cold accretion, indicating 
the Galaxy is indirectly obtaining mass in this fashion. 

There are several properties of the observed HI clouds in the halo that can be used to 
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constrain the models as they are developed further besides the mass limits. As discussed in 
the results section and shown in Fig. 1, most of the observed clouds are found at latitudes 
below 60° and there are several clusters of smaller clouds in specific directions that may 
represent the preferred accretion axes of the Galaxy. The mean negative Vqsr of the cloud 
population (-43 km s"^) is suggestive of a net infall, but the mix of positive and negative 
Vgsr clouds found throughout the sky suggests it is not a simple model. Since we are 
only measuring one component of the cloud's velocity, some clouds with negative or positive 
Vgsr's may be moving away or towards the disk, respectively. In any case, after forming 
the clouds appear to be on a variety of orbits, which after collisions and ram pressure drag 
eventually lead to infall (MB04). There is no correlation between the GSR velocity of the 
cloud and distance if all of the clouds are given similar masses, indicating that the clouds 
do not all have similar masses or reflecting the complex motions of the halo clouds. If the 
ongoing distance determination programs continue to place the HVCs within 60 kpc of the 
Galactic disk, explanations will need to be made for the reason the HI is only found out to 
this radii and what percentage of the baryonic mass the HI represents. Future HI surveys 
being completed by the Galactic Arecibo L-Band Feed Array (GALFA) consortium (e.g., 
Stanimirovic et al. 2006) and the Galactic All-Sky Survey (GASS; McClure- Griffiths et al. 
2006), will be important for characterizing the properties of the clouds and their relationship 
to the Galactic disk. The GALFA surveys, with their increased sensitivity and resolution, 
will be particularly important for examining the flux distribution function and assessing 
how the condensed clouds interact with the diffuse hot halo as they are assimilated into the 
Galactic disk. This analysis will also allow for an estimate of the mass of the elusive diffuse 
hot halo. 

Thanks to James Bullock, Ari Mailer, and Jespcr Sommcr-Larscn for very useful discus- 
sions, to the referee for insightful comments, and to Bart Wakker for providing an updated 
version of the WvW91 catalog. 
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Fig. 1. — Spatial distribution of the 582 potentially condensed high velocity clouds on the sky 
in Galactic coordinates with the negative Vgsr clouds labeled as triangles and the positive 
Vgsr clouds labeled as stars. 
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Fig. 2. — The total mass of the entire population of potentially infalling clouds if they are 
all given the typical distance on the x axis and are assumed to be 10.5% observable neutral 
hydrogen (see text). The HI flux used to calculate this total mass is dominated by some of 
the large complexes like Complex C, so it is more realistic to place the clouds at a distribution 
of distances which keeps the individual cloud masses below 10^ Mq as discussed in the text. 
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Fig. 3. — An example of the total masses of individual clouds if they are assigned a random 
distribution of distances below 150 kpc, are 10.5% observable neutral hydrogen, and are 
confined to have a total mass below 10^ Mq. The total mass in condensed clouds with this 
data is typically 1.1 - 1.3 x 10^ M©. 



